PAIRS OF FREQUENCY-BASED NONHOMOGENEOUS DUAL WAVELET 
FRAMES IN THE DISTRIBUTION SPACE 
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Abstract. In this paper, we study nonhomogeneous wavelet systems which have close relations 
to the fast wavelet transform and homogeneous wavelet systems. We introduce and characterize 
a pair of frequency-based nonhomogeneous dual wavelet frames in the distribution space; the pro- 
posed notion enables us to completely separate the perfect reconstruction property of a wavelet 
system from its stability property in function spaces. The results in this paper lead to a natural ex- 
planation for the oblique extension principle, which has been widely used to construct dual wavelet 
frames from refinable functions, without any a priori condition on the generating wavelet functions 
and refinable functions. A nonhomogeneous wavelet system, which is not necessarily derived from 
refinable functions via a multiresolution analysis, not only has a natural multiresolution-like struc- 
ture that is closely linked to the fast wavelet transform, but also plays a basic role in understanding 
many aspects of wavelet theory. To illustrate the flexibility and generality of the approach in this 
paper, we further extend our results to nonstationary wavelets with real dilation factors and to 
nonstationary wavelet filter banks having the perfect reconstruction property. 
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■ 1. Introduction and Motivations 

In wavelet analysis, we often use translation, dilation, and modulation of functions. For a 
function / : R — > C, throughout the paper we shall use the following notation 

fx±Ax):=\X\ l/2 e- inXx f{\x-k) and f X;k := / A;M = \X\ 1/2 f(X ■ -k), x, A, k, n G R, (1.1) 

O \ where i denotes the imaginary unit. In this paper we shall use d G R\{0} as a dilation factor. In 
applications, d is often taken to be a positive integer greater than one, in particular, the simplest 
case d = 2 is often used. 

Classical wavelets are often defined and studied in the time/space domain with the generating 
Q\ \ wavelet functions belonging to the square integrable function space L 2 (R). For d G R\{0} and for 
a subset ^ of square integrable functions in L 2 (R), linked to discretization of a continuous wavelet 
transform (see [2j [91 [291 EO] ) , the following homogeneous wavelet system 

WS(tf) := {V>di ; * : (1.2) 

is generated by the translation and dilation of the wavelet functions in \l/ and has been extensively 
studied in the function space L 2 (R) in the literature of wavelet analysis. To mention only a few 
references here, see [1] [32]. In this paper, however, we shall see that it is more natural to study 
a nonhomogeneous wavelet system in the frequency domain. It is important to point out here 
that the elements in a set S of this paper are not necessarily distinct and h G S in a summation 
means that h visits every element (with multiplicity) in S once and only once. For example, for 
\l/ = . . . ,ip s }, all the functions ip 1 , . . . ,ip s are not necessarily distinct and ip G \P in (11.21) 
means if) = ip 1 , . . . , ip s . 

Most known classical homogeneous wavelet systems WS(\I/) in the literature are often derived 
from scalar refinable functions (f) or from refinable function vectors <fi = ((f) 1 , . . . ,4> r ) T in L 2 (R) 
([21 [291 SO]). Let us recall the definition of a scalar refinable function here. A function or 
distribution on R is said to be refinable (or d-refinable) if there exists a sequence a = {a(k)}k&z 
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of complex numbers, called the refinement mask or the low-pass filter for the scalar refinable 
function 0, such that 

= \d\J2 a (k)(f>(d--k) (1.3) 

with the above series converging in a proper sense, e.g., in /^(M). Wavelet functions in the 
generating set \I> of a homogeneous wavelet system WS(\I/) are often derived from the refinable 
function by 

if,= |d|£V(Jfe)0(d--A;), V e *, (1.4) 

fcez 

where 6^ = {6^(A;)}fc e z are sequences on Z, called wavelet masks or high-pass filters. For the 
infinite series in ( 11.31) and (11.41) to make sense, one often imposes some decay condition on the 
refinable function and wavelet filters a, b^ so that all the infinite series in (11.31) and (11.41) are well- 
defined in a proper sense. Nevertheless, even for the simplest case of a compactly supported scalar 
refinable function (or distribution) with a finitely supported mask a, the associated refinable 
function with mask a does not always belong to L 2 (M). In fact, it is far from trivial to check 
whether E L 2 (M) in terms of its mask a, see [2"U| |2~T] and references therein for detail. One of 
the motivations of this paper is to study wavelets and framelets without such stringent conditions 
on either the generating wavelet functions 0, if) or their wavelet filters a, for if) E ty. 

For / E -^i(ffi), the Fourier transform used in this paper is defined to be /(£) := f R f(x)e~ lx ^dx, 
and can be naturally extended to square integrable functions and tempered distributions. 
Under certain assumptions, taking Fourier transform on both sides of (11.31) and (11.41) . one can 
easily rewrite (11.31) and (II .4p in the frequency domain as follows: 

0(dO = a(O0(O, o-e. t e R (1.5) 

and 

</>(dO=%)0(O, a.e. £ E R, ip E ^ (1.6) 

provided that all the 27r-periodic (Lebesgue) measurable functions d(£) = J2kez a {k)e~ lk ^ and 

similarly b 1 !' are properly defined. In the following, we shall see that it is often more convenient to 
work with (jl.5p and (II .6p in the frequency domain rather than (II .3p and (11.41) in the time/space 
domain. If there exist positive real numbers r and C such that the 27r-periodic measurable function 
a satisfies |1 — a(£)| ^ C|£| T for almost every £ E [— 7r, it] (this condition is automatically satisfied 
with t = 1 if a is a 27r-periodic trigonometric polynomial with d(0) = 1), for a dilation factor 
d such that |d| > 1, then it is easy to see (also c.f. section 3) that one can define a measurable 
function <p such that 

oo 

<p(0:=na(d- J '0, a-e-£eR. (1.7) 

i=i 

Regardless of whether <£> in (11.71) is a square integrable function or not, <p is a well-defined mea- 
surable function obviously satisfying the frequency-based refinement equation (II. 5p with being 
replaced by ip. The function cp is called the (frequency-based) standard refinable function with 
mask d and dilation d. All the wavelet functions if) in (jl.6p with = are also well-defined measur- 
able functions provided that all 6^ are measurable. This motivates us to study refinable functions 
and wavelets in the frequency domain using (I1.5P and (11.61) so that we can avoid some technical 
issues such as the convergence of the infinite series in (11.31) and (11.41) as well as membership in 
L 2 (M) of the generating refinable function and the generating wavelet functions if> E 

For wavelets derived from refinable functions or refinable function vectors, one of the most 
important key features of wavelets and framelets is its associated fast wavelet transform, which is 
based on the following nonhomogeneous wavelet system: 

WSj($;*) :={(f> d j ;k : k E Z, E $} U {^, k : j > J, k E Z,if> E (1.8) 

where J is an integer, representing the coarsest decomposition (or scale) level of its fast wavelet 
transform. In fact, a one-level fast wavelet transform is just a transform between two sets of 
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wavelet coefficients of a given function represented under two nonhomogeneous wavelet systems at 
two consecutive scale levels. Naturally, for a multi-level wavelet transform, there is a underlying 
sequence of nonhomogeneous wavelet systems at all scale levels, instead of just one single wavelet 
system. For a given WSj ($; at some scale J , we shall see in this paper that via the dilation 
operation it naturally produces a sequence of nonhomogeneous wavelet systems WSj($;^/) for 
all integers J with almost all properties preserved. Consequently, it often suffices to study only 
one nonhomogeneous wavelet system instead of a sequence of them. This desirable property of 
nonhomogeneous wavelet systems is not shared by homogeneous wavelet systems. Furthermore, as 
J — > —oo, the limit of the sequence {WSj($; ^)}j 6 z will naturally lead to a homogeneous wavelet 
system WS(\I/). Hence, in certain sense, a homogeneous wavelet system could be regarded as the 
limit system of a sequence of nonhomogeneous wavelet systems. See section 3 for more detail. 

For a homogeneous wavelet system WS(\I/) that is derived from a refinable function or a refinable 
function vector, due to the absence of a refinable function cf) in the system, the homogeneous wavelet 
system WS(\I/) does not automatically correspond to a fast wavelet transform without ambiguity. 
In fact, the wavelet functions in \1/ of WS(\I / ) could be derived from many other (equivalent) 
refinable functions, which correspond to different fast wavelet transforms with different sets of 
wavelet filters. More precisely, for a 27r-periodic measurable function 9 such that ^ for 

almost every (gM, define := #(£)</>(£), then it is evident that <p is also refinable and satisfies 

km = [0(dOa(O/0(O]fe) and ^(dO = fno/omh), a.e. eei^ef. 

Such a change of generators from a refinable function <p to another equivalent refinable function 
4> is in fact the key idea in the oblique extension principle (OEP) in [31 [TTj, [121 123] (also see 
[131 [HI EH EH 123 [2SI |27] ) to construct compactly supported homogeneous wavelet frames WS(\I/) 
in L 2 (R) with high vanishing moments derived from refinable functions and refinable function 
vectors. See [121 [23] for a detailed discussion on a fast wavelet transform based on a homogeneous 
wavelet system obtained via OEP from refinable function vectors. The effect of the change of 
generators on its fast wavelet transform is addressed in [23]. As we shall see in sections 3 and 
4, nonhomogeneous wavelet systems are closely related to nonstationary wavelets (see [If \7\ I2E]) 
and are naturally employed in a pair of nonhomogeneous dual wavelet frames in a pair of dual 
Sobolev spaces introduced in [27]. Due to these and other considerations, it seems more natural 
and more important for us to study nonhomogeneous wavelet systems WSj($; rather than the 
extensively studied homogeneous wavelet systems WS(\I/). This allows us to understand better 
many aspects of wavelet theory such as a wavelet filter bank induced by OEP and its associated 
wavelets in the function setting without a priori condition on the generating wavelet functions. 

Following the standard notation, we denote by f^(R) the linear space of all compactly supported 
C°° (test) functions with the usual topology, and f^'(R) denotes the linear space of all distributions, 
that is, f^'(R) is the dual space of @(M). By duality, the definition in (11. ip for translation, dilation 
and modulation can be easily generalized from functions to distributions. Moreover, when / in 
(11. ip is a square integrable function or more generally a tempered distribution, for A ^ 0, we have 

f\;k,n = e~ l n f\- 1 ;-n,k an d f\-k = f\-\0,k- (1-9) 

In this paper, we shall use boldface letters to denote functions/distributions (e.g., f, g, ip, if), a, b) 
or sets of functions/distributions (e.g., St) in the frequency domain. 

Let <E» and St be two sets of distributions on R. For an integer J and d E R\{0}, we define a 
frequency- based nonhomogeneous wavelet system FWSj(3>; to be 

FWSj($; *) := {<p d -j ;0 , k : k E Z, if E U {Vd-w : j ^ J, k E Z,ip E i 5f}. (1.10) 

Similarly, a frequency-based homogeneous wavelet system is defined as follows: 

FWS(*) := {VdW : j e Z,fc e Z,<0 e (1.11) 

By f ll.9p . it is straightforward to see that under the Fourier transform, the images of WSj($; ^) 
and WS(^) with $, ^ C L 2 (R) are simply FWSj($; *) and FWS(*), respectively, where $ : = 
{(f) : E $} and : ip E 
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For 1 ^ p ^ oo, by L l ° c (M>) we denote the linear space of all measurable functions / such that 
J K \f\ p <oo for every compact subset K of R, with the usual modification for p = oo saying that 
/ is essentially bounded over K. Note that L l ° c (W) is just the set of all measurable functions that 
can be globally identified as distributions. So, L l ° c (M) is the most natural space for us to study 
wavelets and framelets in the distribution space. For / G L 2 (R), it is evident that / G L 2 (R) C 
L 2 oc (R) C L[ oc (R) C &'(R). However, a distribution if) may not be a tempered distribution 
and therefore, Fourier transform may not be applied so that if) = ip holds for some tempered 
distribution if). Under the setting of tempered distributions on which the Fourier transform can 
apply, although all the definitions and results of this paper in the frequency domain could be 
equivalently translated into the time/space domain by the inverse Fourier transform on tempered 
distributions, to avoid notational complexity and to avoid the a priori underlying assumption that 
/ is a tempered distribution if the notion / is used, it seems very natural and convenient for us 
to work in the frequency domain in this paper. 

For f G ^(R) and if) G Lj° c (R), we shall use the following paring 

(f,if)}:= / mW)di and (</>, f> := (f^) = / if)(OWH- (1.12) 

When f G f^(R) and if) G f^'(R), the duality pairings (f,if)) and (if),f) are understood similarly 
as (f, if>) := (if), f) := if)(f)- Now we are ready to introduce the key notion in this paper. Let 

* = {cp 1 ,...,ip r }, * = {if)\...,if) s } and 9 = {&,..., 9 = {j,\...,fr} (1-13) 

be subsets of @'{R). Let J G Z and d G R\{0}, we say that the pair (FWSj(#; *), FWSj(*; *)), 
where FWSj(3>; ^) is defined in (11.101) . forms a pair of frequency-based nonhomogeneous dual 
wavelet frames in the distribution space fi^(R) if the following identity holds 

r oo s 

E E( f ' Vd-'io,*) (^S-j ; o,*. g) + E E E( f ' ^;o,k) (^d-i ; o,*> g) 

t=i k&L j=j i=i k&L (1.14J 

= 27r(f,g) Vf,gG^(R), 

where the infinite series in (11 . 14j) converge in the following sense: 

(1) For every f , g G ^(R), the following series 

^( f >Vd-^0,fe)(^d-^0,fc>g> aIld 5^^ f '^d-i ;0j ifc)(^d-,';0,fe>g) ( L15 ) 

fcez fcez 

converge absolutely for all integers j ^ J, £ = 1, . . . , r, and £' = 1, . . . , s. 

(2) For every f , g G f^(R), the following limit exists and 

( E E( f ' v5-W (<pUo, k > g) + E E E< f > ^ :n,-) (tfS-^. g>) = Mf , s). 

00 \e=i feez j=j £=i fcez / 

(1.16) 

As we shall discuss in section 3, the above introduced notion enables us to completely separate 
the perfect reconstruction property in (j!.14p from its stability property in function spaces. Since 
the test function space ^(R) is dense in many function spaces, one could extend the perfect 
reconstruction property (or "wavelet expansion" ) in (j!.14p to other function spaces, provided that 
the involved wavelet systems have stability in these function spaces. Let us give a simple example 
here to illustrate this connection for the particular function space L 2 (R). Let <fr and in (11.131) be 
two subsets of distributions in ^'(R). We shall see in section 3 that (FWSj($; *), FWSj($; *)) 
forms a pair of frequency-based nonhomogeneous dual wavelet frames in the distribution space, if 
and only if, <fr, St are subsets of L 2 (R) and 

r oo s 

EEK f ^Uo, fc )| 2 + EEEK f '^-;0, fc )| 2 = 27r ll f llLw Vf GL 2 (R). (1.17) 

i=i fcez j=j i=i 
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Moreover, when |d| > 1, as a direct consequence of (11.171) . one automatically has 

EEE K f '^ ; o,,)| 2 = 2vr||f||i 2(M) V f G L 2 (R). 

jez i=i fcez 

Nonhomogeneous wavelet systems also have a close relation to refinable functions and refinable 
function vectors. Suppose that ^= FWSj(<fr; VP) (that is, multiply every element in FWSj(<fr; \£) 

by the factor -|=) is an orthonormal basis of L 2 (R). Denote (p := (ip 1 , . . . , <p r ) T ■ Without 
assuming in advance that ip is a refinable function vector and all if) 1 , . . . , if) s are derived from 
we can deduce that tp must be a refinable function vector and all if) 1 , . . . , if) s must be derived from 
via similar relations as in (jl.5p and (11.61) . See section 3 for more detail. 

To have some rough ideas about our results on nonhomogeneous wavelet systems, here we 
present two typical results. The following result is a special case of Theorem [61 

Theorem 1. Let d be an integer such that |d| > 1. Let <fr, \I>, <fr, \I> in (11.131) be subsets o/L 2 oc (M). 
Then (FWSj(<fr; SI/), FWSj(3>; \I/)) forms a pair of frequency-based nonhomogeneous dual wavelet 
frames in the distribution space f^'(R) for some integer J (or for all integers J), if and only if, 
the following three statements hold 

(i) For all integers k G Z ; 

r s r 

^^W^(d(e + 27rA;))+^^^^(d(e + 27rA;)) = ^^(0^(^ + 27rA;), a.e. £ G R; (1.18) 
i=i i=i i=i 

(ii) For all integers ko G Z\[dZ], 

r s 

Yy^)^ + 2^ ) + E^W(£ + 2vrA; ) = 0, a.e. £ G R; (1.19) 

(iii) T/ie following identity holds in the sense of distributions: 

r 

lim V^(d-J-)^(d^'-) = 1, (1.20) 



more precisely, lim^+oo ^^ =1 (y? £ (d~' J '-)^(d = (l,f) /or a// f G ^(R). 

In the following, we make some remarks about Theorem [TJ We assumed in Theorem [1] that 
all the generating functions in ^ are from the space L 2 oc (R). Note that L 2 oc (R) includes 

the Fourier transforms of all compactly supported distributions and of all elements in all Sobolev 
spaces. This assumption on membership in L 2 oc (R) can be weakened and is only used to guarantee 
the absolute convergence of the infinite series in (11.151) . See the remark after Lemma [3] in section 2 
for more detail on this natural assumption. 

If we assume additionally that <p e ip e g L 0O (R) for all £ = 1, . . . , r and if (11.201) holds for almost 
every £ G R, by Lebesgue dominated convergence theorem, then f l 1 . 2 f) holds in the sense of 
distributions. If all elements in <fr, <fr, ^ are essentially nonnegative measurable functions, then 
it is not difficult to verify that the conditions in items (i) and (ii) of Theorem [1] are equivalent to 
the following simple conditions: 

r s r 

tfm&m + E tfw&w = E <^ a.e. z e r (1.21) 

1=1 i=i 1=1 

and 



1.22) 



+ 27 rk) = and ^'(0^'(C + 2vrA;) = a.e. £ G R, 
V k G Z\{0}, t = l,...,r, e = l,...,s. 

As we shall see in section 2, items (i) and (ii) of Theorem [1] correspond to a natural multiresolution- 
like structure, which is closely linked to a fast wavelet transform. The condition in item (iii) of 
Theorem [1] is a natural normalization condition which is related to (I1.16p . 
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Comparing with the characterization of a pair of homogeneous dual wavelet frames in the 
space L 2 (R) or a homogeneous orthonormal wavelet basis in L 2 (R) (e.g., see [TUJ HZJ [18j [28j 
I2H [32]), Theorem [1] has several interesting features. Firstly, the characterization in items (i)- 
(iii) of Theorem []] does not involve any infinite series or infinite sums; this is in sharp contrast 
to the homogeneous setting in L 2 (R). Secondly, as we shall see in section 2, all the involved 
infinite sums in the proof of Theorem [1] are in fact finite sums. This allows us to easily generalize 
Theorem [T] to any real dilation factors and to nonstationary wavelets, see sections 2 and 4 for 
detail. Thirdly, we do not require any stability (Bessel) property of the wavelet systems, while 
the homogeneous setting in L 2 (R) needs the stability property to guarantee the convergence of 
the involved infinite series. Fourthly, we do not require in Theorem [T] that the generating wavelet 
functions possess any order of vanishing moments or smoothness, while all the generating wavelet 
functions in the homogeneous setting require at least one vanishing moment. Lastly, from a 
pair of nonhomogeneous dual wavelet frames in L 2 (R), we shall see in section 3 that one can 
always derive an associated pair of homogeneous dual wavelet frames in L 2 (R). In fact, most 
homogeneous wavelet systems in the literature are derived in such a way. We mention that weak 
convergence of wavelet expansions has been characterized in [IB] for homogeneous wavelet systems. 
Similar weak convergence of wavelet expansions that are related to (11.161) also appeared in the 
study of homogeneous dual wavelet frames in L 2 (R) and their frame approximation properties, for 
example, see [T21 CEO US] . We also point out that the approach in this paper can be extended to 
frequency-based homogeneous wavelet systems in the distribution space £^'(R\{0}). 

The following result generalizes the Oblique Extension Principle (OEP) and naturally connects 
a wavelet filter bank with a pair of frequency-based nonhomogeneous dual wavelet frames in & (R) . 

Theorem 2. Let d be an integer such that |d| > 1. Let a, & 1 , . . . , 6 r , b 1 , . . . , b s and a, G 1 , . . . , 6 r , 

b 1 , . . . , b s be 2tt -periodic measurable functions on R. Suppose that there are measurable functions 
(p, (p satisfying 

<p(d£) = a(£)y>(0 and £(d£) = a(£)£(0, a.e. £ G R. (1.23) 
Define * as in fll.!3j) with 

¥/(£):= 0*(£M£), tp e (0--=O e (Om, l=l,...,r (1.24) 

and 

v/(do : = b^(o, := e=i,...,s. a.25) 

Assume that all the elements in belong to L 2 oc (R). Then (FWSj($; *), FWSj($; *)) 

forms a pair of frequency-based nonhomogeneous dual wavelet frames in the distribution space 
f^'(R) for some integer J (or for all integers J), if and only if, 

lim B(d _ - 7 -)(^(d _ J-)(^(d - ''-) = 1 in the sense of distributions (1-26) 

with 

r 

1=1 

and the following fundamental identities are satisfied: 

s 

e(d0^a(0 + ^W)b € (0 = ©(0, a.e.£e<7„n<7£, (1.28) 

e=i 

and 

s 

0(d£MOa(£ + 2-) + W)b'(£ + 2ss) = , a.e. £ G <r v n - 2ss), (1.29) 

e=i 

for alluj = 1, . . . , |d| - I, where - *f - := {£ - ^ : £ G cr^} and 

a v := {£ G R : £ |<p(£ + 2ttA;)| ^ o}, ^:={(GR:^ + 2ttA;)| ^ o}. (1.30) 

fcez fcez 
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In particular, if all 6 l , . . . , 6 r , b 1 , . . . , b s , l , . . . , 6 r , b 1 , . . . , b s are In -periodic measurable junc- 
tions in L l 2 C (WL) and if there exist positive real numbers r and C such that 

|l-a(OKC|e and |1 - a(£)| < C\£\ r , a.e. £ e [-tt, tt], (1.31) 

inen ine frequency- based standard refinable measurable functions ip,<p with masks a, a and the 
dilation factor d , which are defined by 

oo oo 

V (fl:=na(d- J '£) «nd := J[ a(<T'£), (1.32) 

i=l i=i 

are well-defined for almost every £ e R and m /act ip,<p E L^ C (R) . Tnen a// elements in ^ 
belong to L^R). Moreover, (FWSj(«fr; VP), FWSj(<fr; Vt)) forms a pair of frequency-based non- 
homogeneous dual wavelet frames in the distribution space f^'(R) for some integer J (or for all 
integers J), if and only if, the identities ( 11. 281) and ( 11.291) are satisfied for all uj = 1, . . . , |d| — 1, 
and lim^+oo 6 (d~-?-) = 1 in the sense of distributions. 

Note that fll.31f> is automatically satisfied with r = 1 if a and a are 27r-periodic trigono- 
metric polynomials with a(0) = a(0) = 1. A similar result to Theorem [2] also holds when 
ip = (ip 1 , . . . , <p r ) T and <p = (ip 1 , . . . , (p r ) T are refinable measurable function vectors. The identi- 
ties in (11.281) and (11.291) with d = 2 and r = 1 are called the oblique extension principle in |12j . 
provided that all elements in ^ belong to L2(K) and satisfy some technical conditions 

to guarantee the Bessel (stability) property of the homogeneous wavelet systems FWS(^) and 
FWS(*) in the space L 2 (R) (see dH [HI E21 ED E2J). In contrast, our results here generally do 
not require any a priori condition on the generating wavelet functions and provide a natural expla- 
nation for the connection between the perfect reconstruction property induced by OEP in (11.281) 
and (I1.29P in the discrete filter bank setting to wavelets and framelets in the function setting. 

The structure of the paper is as follows. In order to prove Theorems [1] and [21 we shall introduce 
some auxiliary results in section 2. In particular, we shall provide sufficient conditions in section 2 
for the absolute convergence of the infinite series in (|1.15|) . Then we shall prove Theorems [T] and [2] in 
section 2. To explain in more detail about our motivation and importance for studying frequency- 
based nonhomogeneous wavelet systems, we shall discuss in section 3 nonhomogeneous wavelet 
systems in various function spaces such as L 2 (R) and Sobolev spaces, as initiated in [27]. We 
shall see in section 3 that under the stability property, a pair of frequency-based nonhomogeneous 
dual wavelet frames can be naturally extended from the distribution space to a pair of dual 
function spaces. In section 3, we shall also explore the connections between nonhomogeneous and 
homogeneous wavelet systems in the space L 2 (R). To illustrate the flexibility and generality of 
the approach in this paper, we further study nonstationary wavelets which are useful in many 
applications, since nonstationary wavelet filter banks can be implemented in almost the same way 
and efficiency as a traditional fast wavelet transform. However, except a few special cases as 
discussed in (U [7J, [2S] , only few theoretical results on nonstationary wavelets are available in the 
literature, probably partially due to the difficulty in guaranteeing the membership of the associated 
refinable functions in L 2 (R) and in establishing the stability property of the nonstationary wavelet 
systems in L 2 (R). In section 4, we present a complete characterization of a pair of frequency-based 
nonstationary dual wavelet frames in the distribution space. Though the statements and notation 
in section 4 on nonstationary wavelets seem a little bit more complicated comparing with the 
stationary case in sections 1-3, it is worth our effort to provide a better picture to understand 
nonstationary wavelets, since there are few theoretical results on this topic in the literature. 

To understand and study wavelet systems in various function spaces, it is our opinion that there 
are two key fundamental ingredients to be considered. One ingredient is the notion investigated 
in this paper of a pair of frequency-based nonhomogeneous dual wavelet frames in the distribution 
space which enables us to completely separate its perfect reconstruction property from its stability 
property in function spaces. The other ingredient is the stability issue of nonhomogeneous wavelet 
systems in function spaces which we didn't discuss in this paper but shall be addressed elsewhere. 
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2. Frequency-based Nonhomogeneous Wavelet Systems in the Distribution Space 

In this section, we study pairs of frequency-based nonhomogeneous dual wavelet frames in the 
distribution space. To prove Theorems [1] and [2], we first present some sufficient conditions for the 
absolute convergence of the infinite series in fll . 15j) . 

For 1 ^ p ^ oo, by L P (T) we denote the set of all 27r-periodic measurable functions / such that 
\f(x)\ p dx < oo (with the usual modification for p = oo). 

By the following result, we always have the absolute convergence of the infinite series in (11.151) 
provided that all the frequency-based wavelet functions are from the space L 1 ™ 1 



Lemma 3. Let A be a nonzero real number and let ip, ip G L 2 oc (IR). Then for all f , g G @(M), 

J2^,^X;o, k )(^X;o,k,s) = 2tt f £ f (£)g(£ + 2<K\^k) </>(A£)^(A£ + 2itk)dt (2.1) 

with the series on the left-hand side converging absolutely. Note that the infinite sum on the 
right-hand side of (12. ip is in fact finite. 

Proof. By L 0O)C (M) we denote the linear space of all compactly supported measurable functions in 
Loo(R). Note that @(R) C L 00>C (R). More generally, we prove for f , g G L^R). Denote 



h(0:=^f(A- 1 (e + 2 7 rA;))^(e + 27rfc) and h(£) := £ g^^tf + 2nk))^ + 2nk). 

fcez k&L 

Now we show that h, h are well-defined functions in L 2 (T). In fact, since f G L 00jC (R), f has 
compact support and therefore, f is essentially supported inside [— 7TCf , TTCf] for some Cf > with 
Cf depending on f . Now it is easy to see that 

h(0= £ {(\- 1 (Z + 2nk))^ + 27rk), £g[-7T,7t]. (2.2) 

fceZ,|A[^(l+[A|c)/2 



Since ip G L 2 oc (M) and f G Loo, c (R), we see that f(A~ x (- + 2tt£;))V>(- + 2vrA;) G L 2 (M) for every 
k G Z. Therefore, h is a well-defined 27r-periodic function in L 2 (T). Similarly, we have h G £ 2 (T). 
Note that 

(f,Vw) = |A| 1/2 / f (O^Af)e <fc *de = |A|- 1/2 / f (A^O^COe^de = |Ar 1/2 [* h(t)e"*d£ 
Jm. Jr J -it 

and (g, ^A;o,fc) — |A| -1 ^ 2 h(£)e 4fc ^(i£. Since h, h G L 2 (T), by the Parseval identity, we have 
l(f, V'WI 2 = 27TIAI- 1 f |h(0| 2 de < oo, £ |(g, ^ A;0ifc )| 2 = 27TIAI- 1 f |h(£)| 2 ^ < oo, 

and 

with the series on the left-hand side converging absolutely. By the finite sum in ( 12.21) . we have 



iai- 1 / h(0h(e)de=iA|- 1 / f(A- i o^(Ohm= / ^ovkaomao^ 

J2 f(0g(e + 2vrA- 1 fc)^(A0V'(Ae + 27rA;)de, 



fcez 

which completes the proof. □ 

The condition if), if) G L 2 oc (IR) in Lemma [3] is only used in this paper to guarantee the identity 
( 12. ip . As long as (12. ip holds for f , g G f^(R) and the frequency-based wavelet functions belong to 
L °°(B.), all the claims in this paper still hold. Note that L° C (M) is just the set of all measurable 
functions that can be globally identified as distributions. So, L ° C (M) is the most natural and 
weakest space for us to study wavelets and framelets in the distribution space. The condition 
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rjj, ip G L 2 oc (IR) in Lemma [3] could be replaced by other conditions. For example, for every positive 
integer k, if there exist positive numbers Tk and C& such that Tfe > 1/2 and |V>(0 ~~ VKOI ^ 
CfelC - C| Tfc for all f , C G [-&, fc], then it is not difficult to check by ([22]) that h with f G 0(R) is 
a 27r-periodic Lipschitz function with some Lipschitz exponent r > 1/2. By Bernstein Theorem, 
h has an absolutely convergent Fourier series. Now for any if> G L l ° c (M), it is easy to prove that 
( 12. ip indeed holds for all f, g G @(M). Other assumptions could be used to guarantee (12. ip . But 
/^(K) is a large space containing the Fourier transforms of all compactly supported distributions 
and of all elements in all Sobolev spaces. For simplicity of presentation, we shall stick to the space 
L 2 oc (M) for our discussion of frequency-based wavelets and framelets. 

Lemma 4. Let {Xj}°^j be a sequence of nonzero real numbers such that \imj^ +00 Xj = 0. Let 
(p 3 ' 1 , . . . , <p 3,ri and tpi' 1 , . . . , tp 3 '^ be elements in L l £ c (W) with rj G N and j ^ J. Then 



ft EE( f '^)(^oi.g>=2-(f,g) Vf,gG^(R), (2.3) 

i=i fcez 

if and only if, 

r i 

lim y ipi' e (\j-)(p 3,e (\j-) = 1 in the sense of distributions. (2-4) 
i=\ 

Proof. By Lemma [3], we have 

Tj Tj 

EE^<;0, fc )(^ ; o^g)=2vr / £ £ f (0s(£ + ^W)*^(M + 2nk)d£. (2.5) 

£=i fcez ^ R £=i kez 

By f , g G ^(M), f and g are compactly supported. Since lim^+oo Xj = 0, there exists an integer 
J fig such that f (f)g(f + 2ttAt 1 A;) = for all £ G R, fc G Z\{0}, and j ^ J fjg . That is, for 
j ^ max(J, Jf ig ), (I2.5P becomes 

EE( f -4)(4'g) = 27r / f(OgME^W)^(A,o^. (2.6) 

£=i fcez ■' R £=1 

If ( 12 .4p holds in the sense of distributions, then it follows directly from (I2.6P that (12.31) holds. 

Conversely, we can take g G f^(R) such that g takes value one on the support of f, now it 
follows from Q and flUD that 

^Wm ti ^ +oo2?r ti^ 
Hence, (12.41) holds in the sense of distributions. □ 

In the next auxiliary result, we shall study a multiresolution-like structure. More precisely, we 
have the following result. 

Lemma 5. Let X be a nonzero real number. Let ip 1 , . . . , <p r , if} 1 , . . . , if) s , 77 1 , . . . , 77* and tp 1 , . . . , <p r , 
ip 1 , . . . , tp s , r) 1 , ■ ■ ■ ,r) t be elements in L^ili). Then 



£=1 fcez 

EE( f ^W(^;0, fc ,g) Vf,gG^(I 



=1 fcez £=1 fce2 

t 



(2.7) 



=1 fcez 
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if and only if 

4(0 + 4(0 = 4*( A 0, a.e. £ G R, V k G Z n [A _1 Z] , (2.8) 

4(0 + 4(0 = 0, o.e. f G R, V fc G ZVA^Z], (2.9) 

7^(0 = 0, a.e. £ G R, V k G [A^Z^Z, (2.10) 

where A~ X Z := {A _1 fc : k G Z} and /^ fc (0 := £j =1 rfi&tf + 2ttAA;) ; jfc G A _1 Z, and 

4(0:=E^(o^(^ + 27rA; )' 4(0:=E^(o^ + 27rA; )' fceZ - ( 2 - n ) 

1=1 i=\ 
Proof. By Lemma all the infinite series in (12.7ft converge absolutely and (12.71) is equivalent to 

/ E f (Og(£+ 2 ^)(4(o + 4(o)^= / E f(Og(e+2^)4 fe (AO^, 

which can be easily rewritten as 

/ e f (Os(e + 2ttA;) (4(0 + 4(o - /j*(Af))de 

M ^zn[A-iz] 

+ / E f (Og(C + 27rfc) (4(o + 4(o)^e= / E f(Og(e + 27rfc)4 fc (AO^. 

Sufficiency If (12. 8p . (12. 9p . and (12.101) are satisfied, then it is obvious that (12.121) is true and 
therefore, ( 12. 7p holds. 

Necessity. Denote A := ZU [A _1 Z]. For a point x G R, we define dist(x, A) := inf^A \x — y\. By 
(12. 7p . (I2.12p holds. Let k G Z R [A -1 Z] and £o G R be temporarily fixed. Then it is easy to check 
that £ := |dist(&o, A\{fc }) > 0. Consider all f , g G ^(R) such that the support of f is contained 
inside (£o — £, £o + e) and the support of g is contained inside (£o — 2Trk — e, £o — 27r/c + £)• Then 
it is not difficult to verify that 



f (Og(f + 2irk) = o veet,fceA\{M, 

from which we see that ( 12. 12ft becomes 

/ f(Og(e+2^o)(4 o (o + 4°(0 - O a ok = o (2.13) 

for all f , g G ^(R) such that supp f C (( - e, £ + e) and supp g C (( - 2nk — e, £ — 27rfc + £)• 
From (12TT3|) . we must have Jj>(£) + 4°(0 ~ 4^0 = for almost every £ G (£ ~ £,fo + e)- 
Thus, ( 12. 8p must be true. (12.91) and (12.101) can be proved by the same argument. □ 

Now we have the generalized version of Theorem [1] with a general real dilation factor d. 

Theorem 6. Let 6 be a real number such that |d| > 1. Let <fr, <fr, Sf' in (11.131) be subsets of 
L l ° c (R). Then (FWSj($; *),FWSj(l>; *)), where FWSj($; *) is defined in (fTTOj) . /orms a pair 
of frequency-based nonhomogeneous dual wavelet frames in the distribution space f^'(R) for some 
integer J (or for all integers J ), if and only if, 



r 



= 1 

and 



lim \ (p e (d i-){p z (6 = 1 in the sense of distributions (2-14) 



E <^(0^(£ + 2vrfc) + E ^(0^(£ + 2ttA:) 

r 

= E^ I (d 3T 0^(d" 1 (^ + 27rA:)), a.e. £ G R, Vfc € ZD [dZ], 



£=1 



(2.15) 
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r s 

Yy^)^ + 27r *o + £^W(£ + 2nk ^ = °> a - e - ^ R - Vfce z \[ dz ]> ( 2 - 16 ) 

£=i e=i 

r 

£ <^(0^(£ + 2^"^) = 0, a.e.f G R, V jfe G [dZ]\Z. (2.17) 
Proof. By the following simple observation, we have 

(fA;n, fc ,gA;n,fc) = (f,g) and (f A;0 , Q , V>A';n,fe> = ( f , ^A'A^n,*) > A, A' G R\{0}, fc, n G R. (2.18) 

Now it is straightforward to see that for all f , g G 



£ £( f ' ( PU;0,k)( ( PU;0^ g> + £ £( f ' ^d-i;0,fe>(^d-i;O,fc: g) 

r 

£ £^ f ' < ^d-^- 1 ;0 J fe)(^d-i-i;0,fe> g)> 



(2.19) 

£( f > Vd-i-l;0 J fe)(^d-i-l;0,fe>g)> 

£=i fcez 

if and only if, for all f , g G £^(R), 

r s r 

£ £< f > *>i;0, fe ><$;o,*> S) + £ £< f > ^;o jfe )(^;o )fc) g> = £ £<f , rf-W<#--o,*, 6>- (2-20) 

Sufficiency By Lemma with A = d _1 , we see that (I2.20p holds and therefore, (12.191) holds for 
all f , g G 0(R) and j G Z. For J' > J, we define 

r J'-l s 

S J j'(f, g) := £ £< f . VS-W^A*' 6> + E E £< f > ^;o,*> (*-. ; o, fe > g>- (2-21) 
l=\ fcez j=j i=\ kez 

Now by (12.191) . we can easily deduce that 

Si'(f>g) -EE( f ^d--o )fc >(^ i0 ^g)- (2-22) 

i=i kez 

By Lemma HI it follows from (12.141) that limj^ +00 Sj'(i, g) = 2n(f, g) for all f , g G ^(R). Hence, 
(FWSj(3>; \I/), FWSj($; SI/)) forms a pair of frequency-based nonhomogeneous dual wavelet frames 
in the distribution space f^'(R). 

Necessity. By (12.181) . we can easily deduce that (FWSj(<E»; *), FWSj(<fr; *)) forms a pair of 
frequency-based nonhomogeneous dual wavelet frames in the distribution space f^'(R) for some 
integer J if and only if it is true for all integers J. Considering the difference between two 
consecutive integers J and J + 1, we see that (I2.19p must hold and therefore, (12.201) holds. Now 
by Lemma El (j2T5]l . flZD, and (1230]) hold with 77 = 9? and t = r, or equivalently, f[2TT5]) . fl2~T6|) . 
and (I2.17P hold. Since (I2.19P holds, we deduce that (I2.22p holds. By Lemma HI it follows from our 
assumption \im.j/_, +00 Sj (f , g) = 27r(f , g) that (I2.14p must hold. □ 

We point out that Theorem [6] and the approach in this paper can be extended to frequency- 
based homogeneous wavelet systems in the distribution space ^'(R\{0}), which is the dual space 
of the test function space i^(R\{0}) consisting of all compactly supported C°° functions whose 
supports are contained inside R\{0}. We shall address this issue elsewhere. 

Now we are ready to prove Theorems [1] and [21 

Proof of Theorem [IJ By Theorem [H] with d being an integer, it suffices to show that (jl.lBp and 
([USD are equivalent to the three conditions fl2TT5]l . fl2~7TB]) . and fLTT7|) . Note that [dZ]\Z is the 
empty set. So, (12.171) with d being an integer is automatically true. It is evident that (12.161) is 
equivalent to (11.191) . Since d is an integer, we have Z D [dZ] = dZ. Now it is also easy to see that 
(I1.18P is equivalent to (12.151) . This completes the proof. □ 

We use Theorem [1] to prove Theorem [2] as follows: 
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Proof of Theorem^ We prove the first part of Theorem [2] first. By (11.231) and (11.241) . we have 

r r 

i=i e=i 
and similarly by (11.25j) . 

s s 

i=i e=i 
for all integers k G Z. Now fll.l8f> is equivalent to 

s 

^I)y>(£ + 2nk) (Q(dOW)m + E W)b"(0) = ^0£(£ + 27rA;)e(0 (2.23) 

i=\ 

for all fceZ. Now it is not difficult to deduce that (I2.23f) is equivalent to (11.281) . 

Note that any k G Z\[dZ] can be uniquely written as k = to + dk for G {1, . . . , |d| — 1} and 
G Z. Replacing £ in f 1 1 . 1 9 j) by d£, by the same argument, we see that fll . 19j) is equivalent to 

s 

+ 2 -ir + ^k) (e(do^)a(e + 2 -T) + Yl + ¥)) = ( 2 - 24 ) 

for all fc G Z and cj = l,...,|d| — 1. Now it is not difficult to deduce that (12.241) is equivalent to 
( TTT29|) withw = l,...,|d| -1. 

We now prove the second part of Theorem [21 Note that (11.311) implies that a, a G L^R) and 
max(|l-a(d- J 'OUl-a(d-' J 'OI) < C|d| _TJ '|f | T for all ^ G (-|d|%, |d|%]. It is a standard argument 
to show that both ip and <p in (I1.23P are well-defined measurable functions in L^ C (R) and 

lim cp(d- j £) = 1 = lim (p(d~ j O, a.e.^GR. (2.25) 

In fact, by the same argument as in [221 Page 93] or [26, Page 932], (11.311) also implies that for 
any e > 0, there exists c > such that 

l-e < \<p{£)\ < 1+e and 1 - e < \ip(£)\ < 1 + e, a.e. £ G [-c,c]. (2.26) 

See the proof of Theorem [16] in section 4 for more detail on proving (I2.25P and (I2.26p . Take e = 1/2 
in (I2.26p . By the definition of in (11. 27ft . we conclude that 

|9(d-^)^(d^e)^(d-^)| < ||e(d~^)| < 9|6(d-^)Rd^)^(d-^)| (2-27) 

for almost every £ G [— c, c] and for all j ^ 0. Consequently, by the generalized Lebesgue dominated 
convergence theorem and using (12.251) and (I2.27p . we can conclude that (I1.26P holds in the sense 
of distributions if and only if lim^+oo 0(d _:? '-) = 1 holds in the sense of distributions. 

Since <p, <p G L l £(R) and 1 , . . . , 8 r , b 1 , . . . , b s , 1 , . . . , r , b 1 , . . . , b s G L l ° c (R), by the definition 
in (11.241) and (11.251) . it is evident that all the measurable functions in <fr, <fr, ^ belong to the 
desired space L 2 oc (R). By what has been proved for the first part of Theorem [21 the claim in the 
second part of Theorem [2] holds. □ 

3. NONHOMOGENEOUS WAVELET SYSTEMS IN FUNCTION SPACES AND L 2 (R) 

In this section, we shall discuss nonhomogeneous wavelet systems in the space L 2 (R) and other 
function spaces. We shall see that a pair of frequency-based nonhomogeneous dual wavelet frames 
in the distribution space plays a basic role in the study of wavelets and framelets in various function 
spaces. We shall also discuss the connections between nonhomogeneous and homogeneous wavelet 
systems in L 2 (R). 

Let us recall some necessary definitions first. For r 6 K, we denote by if T (R) the Sobolev space 
consisting of all tempered distributions / such that 



ht (r) :=±- [ |/(0| 2 (l + |ei 2 )X<oo. (3.1) 
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Note that H~ T (M.) is the dual space of if r (IR), since / G H~ T (M.) can be regarded as a continuous 



linear functional on H T (R) in the sense of (/, g) — j- L f° r 9 e ^ 

Let d be a nonzero real number and WSj(<3>; be a nonhomogeneous wavelet system defined 
in (11.81) . We define a normed sequence space b_ffr( R ), indexed by the elements of WSj($; \&), with 
weighted norm as follows: 

oo 

||M, e ws.H*;dlL (R) :=£^ (3-2) 

fcez j=J Ve* fcez 

For $, \l/ C H T (W), we say that WSj($; has stability in H T (¥L) with respect to the normed 
sequence space b^r( R ) if there exist positive constants C\ and C 2 such that 

Ci\\f\\ 2 H -T [R) < \\{(f,h)} heWSA *m\\l H T (1l) <C 2 \\f\\ 2 H . T{m V/efT-(R), (3.3) 
or more explicitly, 

oo 

ciimi 2 H -. (R) ^ EE idr 2Tj K/,0d^)i 2 + E E E idr 2 m^ ;fc >i 2 < ^ii/ii^ (r) . (3.4) 

Note that if r (R) is a Hilbert space under the inner product 

(/,ffW) := ^ / /(W + l£lT^ /, <? G H^(R). 
An Jr 

It was shown in [27] Proposition 2.1] that (13.41) holds for all / G H~ T (W), if and only if, 



(3.5) 



Cl|bll|r (R) ^ EEl^'^l TJ M^(«)| 2 + EEEl^l d l T ^;fc)^(R)| 2 

In other words, that WSj($; \l/) has stability in iJ T (M) with respect to 6h t (k) is equivalent to 
saying that after a proper renormalization of WSj($; \l/), 

WS}($;^) :={|d|- Tj d , ;fc : G Z, G $} U {|d|-^. fe :^J,fceZ,^e*} (3.6) 

is a frame in the Hilbert space if T (IR). For ^ G if T (M) and |d| > 1, in the following we show that 

^iiMI^( R ) < ll|dr Ti ^ ;fc ||^( R) < M\i>\\m W V j >J,keZ, (3.7) 

for some positive constants A\ and A 2 depending on d, r, J and ip- F° r t < 0, to prove (13. 7L we 
further assume that 

/ "|^(0rie| 2 X<oo, (3.8) 



which is also a necessary condition for ( 13. 71) to hold. We now prove (13.71) . In fact, 

1 f - /ldl~ 2j ' 4- I£I 2 \t 

iiidr T ^ ;fc iiW) = ^ y R i^)i 2 (i+iei 2 r( ' ) (3.9) 

Since |d| > 1, it is easy to deduce that 

n 2 ^idi^+ier^ui-w 



i7I^<|d| +1, Vj^ J,eeR. (3.10) 



1 + lel 2 

For r ^ 0, noting that there exists e > such that 

^ / i^(0i a (i + ier) T de>5ii^i^ w , 

we deduce from (E2D and (13TTUD that (J3T7]) holds with A x = (1 + £~ 2 )~ T /2 and A 2 = (|d|- 2J + 1) T . 
For r < 0, by our assumption in (13. 8p . there exists e > such that 

^ / l^(OI 2 iei 2r ^^^ll^llW). 

2tt ./^ c ro : | ? |^ e } 2 
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Now we can easily deduce from ( 13. 9ft and ( 13. lOf) that ( 13. 7ft holds with Ai = (|d|~ 2J + 1) T and 
A2 = \ + (1 + e~ 2 )~ T . Using (13. 9p and applying Fatou's lemma to (13. 7p with j — > +00, it is easy 
to see that (13.81) is a necessary condition for ( 13.71) to hold. 

Hence, the weights \d\~ T ^ in the definition of || ■ ||b wr(R) is chosen to be the normalization factors 
of the elements ipdj-k m the space H T (R). The stability property in (13.41) . which is equivalent 
to the frame property in (13.51) . characterizes the Sobolev space H~ T (R). We point out that the 
characterization of homogeneous Besov spaces by homogeneous framelets in [TJ [15] uses nonlinear 
approximation and is quite different in nature to ( 13.41) and ( 13.51) . 

Now we are ready to recall the definition in [27] on a pair of nonhomogeneous dual wavelet 
frames in a pair of dual Sobolev spaces (H T (R), H ~ r (M)). Let d be the dilation factor. Let 

$ = {0\...,0 r }, * = {V>V..,V'} and $ = {0\...,^}, * = {$\. ..,$'} (3.11) 

be subsets of tempered distributions. We say that the pair (WSj($; \I/), WSj(<3>; \I/)), or more 
precisely (WS}($; WSj T ($; \I/)), forms a pair of nonhomogeneous dual wavelet frames in a 
pair of dual Sobolev spaces (H T (R), H~ T (R)) if 

(1) C H T (R) and WSj($; *) has stability in H T (R) with respect to b m{K) . That is, 
WS}($; *) is a frame in H T (R); 

(2) $, * C #~ T (M) and WSj($; #) has stability in #~ r (M) with respect to b ff -r (K) . That is, 
WS7 T ($; *) is a frame in H- T (R); 

(3) for all / G #~ r (M) and g G # T (M), the following identity holds 

r oo s 

<f> 9) = E E</ > + E E E</> ^ ;fe > (3.12) 

1=1 k&L j=J e=i feez 

with the series on the right-hand side converging absolutely. 

The above definition is introduced in [27]. When r = 0, since I/2(M) — H°(R), it is easy to 
see that the above definition of a pair of nonhomogeneous dual wavelet frames in (/^(IR), L2(R)) 
becomes the definition of a pair of nonhomogeneous dual wavelet frames in /^(M). 

Suppose that the pair (WS j($; \1/), WSj($; forms a pair of nonhomogeneous dual wavelet 
frames in (H T (R), H~ T (R)). Using (13.121) and the upper stability (that is, the right-hand inequality 
in ( 13.41) holds) of the two nonhomogeneous wavelet systems, it is not difficult to see that we have 
the following representations: 

r oo s 

f = E E</> <A>M>* + EEEtt ; G ^ T(M) 

1=1 k&Z j=J 1=1 keZ 

with the series converging unconditionally in the space H~ T (R), and similarly, 

r 00 s 

9= E E& %>M<* + E E E ^ * ; *>v& ;fc , <? e # r w 

£=i fcez j=J l=l kez 

with the series converging unconditionally in the space H T (R). See [27] for more detail. 

Let 23 be a normed function space and 23' be its dual. For example, for Besov spaces BZ" (R), 
(-B£,(R))' = where r £ R,l< p < 00 and H g < 00 with l/p + 1/p' = 1/q + 1/q' = 1. 

Similarly, for Triebel-Lizorkin spaces F^ q (R), {F^ q (R))' = F^ ql (R). Replacing H T (R) and H~ T (R) 
by 03 and 03', respectively, the notion of a pair of nonhomogeneous dual wavelet frames can be 
generalized from a pair of dual Sobolev spaces (H T (R), H~ T (R)) to a pair of dual function spaces 
(23, 23') by a proper choice of some normed sequence spaces brg and b<s>. We shall not further 
address this issue in this paper. 

In the rest of this section, we shall use the following notation 

lf(R) := {/ : / G iT(R)} and WfWgz* := / G H T (R). (3.13) 
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By the following result, we see that the notion of a frequency-based nonhomogeneous dual 
wavelet frames in the distribution space plays a basic role in the study of pairs of nonhomogeneous 
dual wavelet frames in a pair of dual function spaces. 

Theorem 7. Let r G R and d be a nonzero real number. Let <£>, $, \P in ( 13. lip be subsets of 
tempered distributions. Define 

$ := {(f) : 0G$}, : ip £ $ := {</> : G $}, * := : $ G (3.14) 

T/ien (WSj($; WSj($; \l/)) forms a pair of nonhomogeneous dual wavelet frames in the pair 
of dual Sobolev spaces (if T (R), H~ T (R)), if and only if, 

(i) there exists a positive constant C such that 



and 



{<f,h>} heF ws. j( *;*)|| fa _ R) < C||f|l^)> f G 0(R) (3.15) 



|{(f>h)k eF ws j( ^)IL^ T(R) < C W f \\my f e ( 3 - 16 ) 



(ii) the pair (FWSj($; *), FWSj($; *)), w/w'c/i zs tfie ma^e of the pair (WSj($; WSj($; #)) 
under the Fourier transform, forms a pair of frequency-based nonhomogeneous dual wavelet 
frames in the distribution space £^'(R). 

Proof. For / G #(R) := {ft:/iG 0(R)} and # G F a (R) with a G R, the relation (/,#) = ^(/,#) 

holds. Since #(R) is contained in both # T (R) and #~ T (R) (or equivalently, 9{R) C iF(R) and 

@(M) C H~ T (§t)), the necessity part is evident. Hence, it suffices to prove the sufficiency part. 
By (BZE5D , for all f G ^(R) and iVeN, 

N N N 

E E idr 2 ^i(f^ d -.w)i 2 + EE E idr^Kf^d-^i^^iifii^. 0.17) 

i^G* k=—N j=J Ve* k=—N 

Using (l3~l7D . we now prove that * C ^(R). By (l3~i7l) . for tjj G we have | (f , Vd-J ; o,o) I < 
C|d| rJ ||f lljj-rw^x for all f G f^(R). Therefore, (•, i/ J d- J ; o,o) can De extended into a continuous linear 

functional on if~ r (R). Since if T (R) is the dual space of i/~ T (R), there exists G iJ r (R) such 
that (•, iftd- J ;0,o) — V*)- m particular, (f , t/>d-J ;0 ,o — V>) = for all f G 5?(R). Since ^d~ J ;0,o ~~ "0 
is a distribution, we must have tpd- J -,oo = "0 in the sense of distributions. By G H T (M), we 
deduce that i^d- J -,o,o -f^ r (R) and therefor, if) G if T (R). Consequently, by (I3.17p . we proved that 
* C H T (R). In other words, we proved that $, * C iJ r (R). 

Since ^(R) is dense in _£/~~ r (R) and since all (■, ip d -j. o k ) , cp G $ and (•, V'd-jjo.fe)? V' £ ^ are 
continuous linear functionals on iJ~ T (R), we see that (13.171) holds for all f G H~ T (M.) and all iV G N. 
Letting iV -> +oo in (13717]) and noting (/,#) = for all / G #~ r (R) and 5 G # r (R), we 

conclude that $, * C IP(R) and 

r 00 s „2 

idr^-^K/-, ^^ ;A; >i 2 -i- 5E idr 2 ^i</-, ^ ;fc >i 2 ^ ^rii/'ii^-^dR), /e#- r (R). ( 3 .is) 

1=1 fcez i=j £=1 k&L 

Similarly, we can show that (13.161) implies $, C H~ T (M.) and 

r 00 s „2 

5^5^ |d| 2 -^|<^,^. ;fc >| 2 + |d|— |<^,^ ;fc >| 2 ^ I^II^II^cm), geir{R). (3.19) 

£=i fcez j=J t=\ kez 

Define two operators W : H~ T (M) — > b H ^(M) and W : if r (R) — > bfl- T (R) by 

W := {(/,/i)}hews. 7 (* ; *) and W# := {(g, tyhewSj&vy ( 3 - 20 ) 
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Then (13. 18ft and ( 13.191) are equivalent to saying that the operators W and W are well-defined 
bounded linear operators, more precisely, 

\\Wf\\ bHTm < -^L\\f\\ H - TW , \\Wg\\ b ^ -^lbll^(K) V/ G H-T(R),ge H T (R). (3.21) 

V Z7T V Z7T 



By Cauchy-Schwarz inequality, using (13. 18j) and (I3.19p . we see that for / G if ~ T (IR) and g G IF 
the series on the right-hand side of ( 13.121) converges absolutely. Now by assumption in item (ii), 
we see that (137L2D holds for all f,g G @(R). In other words, for f,g e 

(Wf,Wg) = {f,g). (3.22) 

Note that ^(M) is dense in both if~ T (M) and H T (R). Now we use a standard argument to show 
that (I3T22|) holds for all / G #" T (M) and # G # r (R). For / G H- T (R) and # G # r (R), there exist 
two sequences {/n}^ and {<?n}5?Li in ^(R) such that 

lim \\f n - f\\ H -r (k) = and lim ||# n - pHh'-(k) = 0. (3.23) 

n— >oo n— ioo 

Observe that 

(W/,Vto) = (W(/ - f n ),Wg) + (Wf n ,W(g - g n )) + (Wf n , Wg n ). 
By f n ,g n G #(R) and (|£22), we have (Wf n ,Wg n ) = (f n ,g n )- Therefore, we have 

W, Wg) - (/, = (W(/ - f n ),Wg) + (W/ n , 1% - g n )) + (/„ - /, g n ) + (/, ^ - g). 
Now by (13.211) and the triangle inequality, 

I (W/, W^) - </, s> I < I (W(f - f n ) , Wg) | + | (W/ n , W(g - g n ))\ + \(f n - /, g n ) \ + | (/, g n - g) \ 

^ Cl (jl/ — /n||if- T (R)||<?||i^(R) + ||/n||ii"- r (R)||5 — 9n 1 1 H T (R) 
+ ||/n - /||tf-'-(R)||fl , n||/r r (R) + 1 1 / 1 1 (R) Hs'n ~ #11^(1 



where Ci := ^- + 1. By (13.231) . the right-hand side of the above inequality goes to as n — > oo. 
Consequently, " k^D holds for all / G #~ r (M) and # G # r (M). That is, (13321 and item (3) have 
been verified. 

To prove item (1), by f!3.22|) and (13.211) . we have 

\(f,g)\ = \(Wf,Wg)\ < ||Wlk n J%lk_. (B) < \\g\\ HT(m ||W/|| bllT(R) . 

That is, for / G H~ T (R), the following inequality holds: 



ll W /l|b^ (K) ^ SUP jr-r = — -||/|| ffT(R)j 

^ ff e/f T (M)\{0} ||<7||i? T (R) <^ 

where in the last step we used the fact that H T (R) is the dual space of H~ T (R). Hence, item (1) 
holds. Item (2) can be proved similarly. This completes the proof of the sufficiency part. □ 

Next, we discuss connections between nonhomogeneous and homogeneous wavelet systems in 
the particular function space L2(R). To do so, we need an auxiliary result, which is essentially 
known in the literature, e.g., see [231 Page 28]. For completeness, we present a proof here. 

Lemma 8. Let d be a real number such that |d| > 1. Let <fi G L20R.) such that there exists a 
positive constant C such that J2kez \{fi4>{' ~ ^))| 2 ^ C|I/IIl 2 (m) f or a ^ f e ^(R) ( or equivalently , 
J2kez \<t>(£ + 27rA; )| 2 < C f or alm ost every £ G R). Then 

lim V|(/,0 dJ;fc )| 2 = O V/gL 2 (R). (3.24) 
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Proof. It suffices to prove the case d > 1, since the negative case d < — 1 can be proved similarly. 
So, we assume d > 1. We first show that (I3.24p holds for all / = X{h,t 2 ) with t\ < t 2 , where X(ti,t 2 ) 
denotes the characteristic function of the open interval (t±, £ 2 ). For / = X(u,t 2 ), by Cauchy-Schwarz 
inequality, we have 

2 



\(f,M\' 



4>(x)dx 



^d'H / dx)( / \(f){x)\ 2 dx). 

' (dH 1 -k,dH 2 -k) ' ^ J {dH\—k,dH 2 —k) 



Since f( d3tl _ k dn 2 -k) ^ x = ^{tz — t±), noting that lim^-oo d J = by d > 1 and all (dHi — k, dH 2 — 
k),k G Z are disjoint as j — > — oo, we deduce that 

^K/,0d,; fc )| 2 ^|t2-tl| / \<P(x)\ 2 dx - 0, 

i.^n? JVu c w(dHi—k,dno—k) 



as j — > — oo, since G L 2 (R). Consequently, (I3.24p holds for all / = X(ti,t 2 )- So, Q3.24H holds for 
all / that are finite linear combinations of characteristic functions of bounded open intervals. 

On the other hand, define operators Pj : L 2 (R) — > Z 2 (Z) by Pjf := {(/, 4>di-,k)}k&- By Pjf = 
Pofd-j;o,o, we have 

H-fj/ll/ 2 2 (Z;) = ll^o/d-J;0,o||fe(Z) ^ C|l/d-J;0,olli 2 (]R) = ^ll / lli 2 (R) • 

When Ylkez l^(£ + 27r A:) | 2 ^ C for almost every £ G R, by Lemma [3l we see that 

2 



ll^/lfiw = ^ E I (/' fe) I' = £ T | E /V'(£ + 2^))0(£ + 2vrA;) 

^r- T (Ei^'^+ 27rA; ))i 2 )(Ei^+ 27rA; )i 2 )^ 
^f dj rEiA di (e+2 7 rA ; ))i 2 rfe = ^ / i/(e)i 2 = cn/f i2(R) . 

Let / G L 2 (R). For an arbitrary e > 0, there is g, which is a finite linear combination of charac- 
teristic functions of bounded open intervals, such that II/— 5 , ||l 2 (R) ^ e - Since lirn 7 _ > _ 00 ||-fjfl , ||i 2 (z) = 
0, there exists J such that ||-P 7 -fl'||i 2 (z) ^ £ for all j ^ J. Hence, for all j ^ J, we have 

WPjfLm < WW - g)\\h<B + ll^llbpo <y/c\\f- ghm + H^IU) < (Vc + i)e. 

Hence, linij^-oo ||-Pj/||; 2 (z) = 0. That is, (I3.24p holds. □ 

Theorem 9. Let d be a real number. Let $, $, \I/ in A3. 1 1 j) 6e subsets o/L 2 (R). Suppose that 
WSj($; which is defined in ( ll.8p . is a frame in L 2 (R) /or some integer J, that is, there are 
positive constants C\ and C 2 such that 

r oo s 

^ill/lli 2W ^EEl^^)| 2 + EEEl^^;^| 2 ^^||/||i 2W , / G L 2 (R). (3.25) 

l=\ k&L j=J 1=1 k&L 

Then (13.251) holds for all integers J. If in addition \d\ > 1, then WS( 1 i r ) must be a frame in L 2 (R) 
with the same frame bounds satisfying 

Ci||/IIL( R ) < EEEl^^>)| 2 < ^||/||i 2(R) , / G L 2 (R). (3.26) 

jez i=\ k&L 

If \d\ > 1 and the pair (WSj($; \l/), WSj($, iff)) forms a pair of nonhomogeneous dual wavelet 
frames in (L 2 (R), L 2 (R)) for some integer J, then (WS(\I/), WS(\I / )) forms a pair of homogeneous 
dual wavelet frames in L 2 (R), that is, both WS(\I/) and WS(\&) are frames in L 2 (R), and the 
following identity holds 

s 

(f>9) = £££</.V&; fc )<#i; fc ,0>, f,9 G ^(R) (3.27) 
jez t=i fcez 
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with the series converging absolutely. 

Proof. By the simple observation in (12.181) . it is easy to see that for all integers J, (13. 25ft holds 
with the same constants C\ and C 2 . In particular, for all J £ Z, we have 



IWII 2 :=££l</^>l 2 ^ 



!l 2 ( 

=i fcez 



Since |d| > 1, by Lemma [BJ we have limj^.oo Y^e=i T^k^i I (/> ^d J -fc) I 2 = ^ ^ or an / e ^(R)- Now 
it is easy to deduce that ( 13.261) holds. 

To prove the second claim, by what has been proved, both WS(\I/) and WS(\I/) are frames in 
L 2 (R). Note that the pair (WSj($; WSj($, \l/)) forms a pair of nonhomogeneous dual wavelet 
frames in (L 2 (R), L 2 (R)) for all integers J. Thus, for all integers J, we have 

r oo s 

V>9) = ££tf + £££(/. ^ ifc >(^ ifc ^>, f,9 G (3-28) 

with the series on the right-hand side converging absolutely. By Lemma [8] again, for f,g £ L 2 (R), 
we have lim J ^_ 00 J2kez(f' ^ ; fc) ($F ; fc> = °- Now we see tnat fl3 - 27 P holds. □ 
As a direct consequence of all the above results, we have 

Corollary 10. Let d be a real number such that |d| > 1. Let $ = {y? 1 , . . . , ip r } and = 
{i/? 1 , . . . , if) s } be subsets of distributions in f^'(R). Tnen ine following statements are equivalent 

(1) FWSj(<fr; VP), which is defined in (ll.lOp . is a frequency-based nonhomogeneous tight wavelet 
frame in L 2 (R) for some integer J, that is, ^ C L 2 (R) and 

r oo s 

£ £ I ( f . <fUo,k)\ 2 + £ £ £ I ( f > *Pd-i;0,k)\ 2 = 27r||f ||i 2(K) V f £ L 2 (R); (3.29) 

(2) item (1) is true for all integers J; 

(3) (FWSj(<fr; \&), FWSj($; forms a pair of frequency-based nonhomogeneous dual wavelet 
frames in the distribution space f^'(R); 

(4) * C L 2 oc (M) and (I2~TD . (I2TT5D . (12TB . (12TT7D hold with $ := $ and * := 

(5) there exist 1 , . . . , r , . . . , -0 s £ L 2 (K) sitc/j taai cp 1 = 1 , . . . , <£> r = </> r , i/? 1 = . . . , i/> s = 
■0 s , and WSj({0 1 , . . . , r }; . . . , ^ s }) a nonhomogeneous tight wavelet frame in L 2 ( 



£ £ I (f, <f>Uk)\ 2 + £ £ £ I (/» v4;*> I 2 = v / £ l 2 (r). (3.30) 

1=1 k& j=J 1=1 kez 

Moreover, any of the above statements implies that FWS(^) is a frequency-based homogeneous 
tight wavelet frame in L 2 (R) ; that is, C L 2 (R) and 

s 

£££i<f,t/4 i0 , fe >i 2 = 2vr||f||i 2(M) Vf £ L 2 (R). (3.31) 

jez £=i fcez 

Proof. We first show that if item (3) holds, then (13.151) must be true with C = In and r = 0. 
Since item (3) holds, by definition, for all f £ £^(R), we have 

,5$. (EEK'.^wi'+EEEK'.^wi') = 

l=\ fcez j=j £=i fcez 

Now it is straightforward to see that (13.151) must be true with C = 2tt and r = 0, since ||/ = 

2^ ll^ll 2 2 (M)- Consequently, by Theorem^, item (3) implies <E», ^ C L 2 (R). 

Now by Theorem [7J (1), (2), and (3) are equivalent to each other. The equivalence of (3) and 
(4) is guaranteed by Theorem [6j The equivalence between (1) and (5) is trivial. (13.311) is a direct 
consequence of Theorem [91 □ 
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From the following result, we see that there is a natural connection between refmable function 
vectors and frequency-based nonhomogeneous orthonormal wavelet bases in L20R). 

Proposition 11. Let A be a nonzero real number. Let $ and ^ in (11.131) be subsets o/L 2 (IR). 
Suppose that FWS *) is an orthonormal basis of L 2 (M.) for some integer J (this is equiva- 
lent to saying that (I3.29P holds with H^Hz^R) = ■ ■ • = ||yj r ||i 2 (R) = ||V' 1 ||l 2 (k) = • • • = ||V1|l 2 (ir) = 
\[2t\). Denote cp := (cp 1 , . . . , <p r ) T and i\) := (if) 1 , ■ . . , ift s ) T ■ Then there must exist r x r matrix a 
and s x r matrix h of 2n -periodic measurable functions in L 2 (T) such that 

£(dO = a(0£(0 and t?(d£) = b(£)£(£), a.e. (el. (3.32) 
Moreover, if \d\ > 1, then ^= FWS^) is also an orthonormal basis 0/L2OR). 

Proof. Note that ^7= FWSj(<fr; SI/) is an orthonormal basis of L 2 (M) if and only if it is an or- 
thonormal basis of L 2 (M) for all integers J. Consider the expansion of the elements <f^. 00 
and V'd-o.o under the orthonormal basis FWSo(<J?; By orthogonality, we have <£>d' ;0 = 

i E<=i £k 6 z<Vd;o,o: Vi;o,fe)Vw and ^d;o,o = ^ELiEitezW^O'^W^W Noting that 
y>f ;0)O (O = \d\ l l 2 <f? {d£) and <p{. 0>k {£) = e- ifc V(0, we see that holds with 

[ a(£) W = ^T^T E^Jo.o, ^i ; o,,)e-^, * = 1, . . . , r 



2^ 

and 



[b(£)W = 7 "TFjf E«o, 0) ¥>i ; o,*>^, * = 1, • • • , r, e = 1, . . • , s, 
27r Vl d l fce z 

where [a(£)]^ denotes the (£', £)-entry of the matrix a(£). Since ^ fcgZ Kv'd-oc 

[a]^ are well defined elements in L 2 (T). Similarly, all \b]t,£ are we U defined elements in L 2 (T). □ 



4. NONSTATIONARY DUAL WAVELET FRAMES IN THE DISTRIBUTION SPACE 

All the results in the previous sections have been mainly built on the multiresolution-like struc- 
ture in (I2.19P for stationary nonhomogeneous wavelet systems. Here stationary means that at scale 
level j the dilation is d J and the generating wavelet functions are independent of the scale level j. 
The result in Lemma [5] characterizing the multiresolution-like structure in (12 .7p makes most proofs 
in the previous sections relatively simple. Nonhomogeneous wavelet systems are closely related to 
nonstationary wavelets, which are useful in many applications since the nonstationary wavelet fil- 
ter banks can be implemented in almost the same way and efficiency as a traditional fast wavelet 
transform. However, except a few special cases as discussed in [U [7J, EE] and some references 
therein, only few theoretical results on nonstationary wavelets are available in the literature. 

In this section, we shall see that the notion of a pair of frequency-based nonhomogeneous dual 
wavelet frames in the distribution space is very flexible and similar results hold in the most general 
setting of fully nonstationary wavelets. Since there are few theoretical results on nonstationary 
wavelets in the literature, it is worth our effort to provide a better picture to understand them in 
this section. 

Let us first introduce the notion of a pair of frequency-based nonstationary dual wavelet frames 
in the distribution space. Let J G Z and {Aj}^j be a sequence of nonzero real numbers. Let 

* = {¥>\. ..,¥>'}, * = {<?,... ,0*} (4.1) 

and 

& ■= {tpi>\ . . . , ip**}, & : = . . . , (4.2) 

be subsets of distributions in f^'(R) with j ^ J and Sj G N. We say that the pair 

(FWSj(#; {&}? =J ), FWSj(#; {&}T=j)) (4-3) 
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forms a pair of frequency-based nonstationary dual wavelet frames in the distribution space &{} 
if the following identity holds: 



OO Sj 



E E< f ' *>U*> §) + E E E< f ' <;0, fc ) (#£0,*. §) = 2*<f , g> V f , g G ^(K), (4.4) 

£=i fcez j=J £=i fcez 

where the infinite series in (14.41) converge in the following sense: 

(1) For every f , g G 0(R), all the series E fce z( f > Va j; o,*> (<£i j; o,fc; g> and E fce z( f > ^o,fc> g) 
converge absolutely for every integer j ' ^ J, I — 1, . . . , r, and £' = 1, . . . , sf, 

(2) for every f , g G @(M), the following limit exists and 



jPZ. (EE( f '^ft Wl g>/EEE( f ^y*^.6> ) =27r(f,g). (4.5) 



The stationary nonhomogeneous wavelet systems considered in previous sections correspond to 
the case that Xj = d~ J , Sj = s, and Vf^ = ^ for all j ^ J; that is, the generating wavelet functions 
remain stationary (unchanged) at all the scale levels j. 

A pair of frequency-based nonstationary dual wavelet frames in the distribution space @'(M) 
can be similarly characterized by the following result. 

Theorem 12. Let J be an integer and {Xj}°^j be a sequence of nonzero real numbers such that 
lim^+cx: = 0. Let <fr, <fr in P~T]) and \I> J , & in (14721) fre subsets of L l £ c (M) for all integers j ^ J. 
JTien t/ie pair in (14.31) forms a pair of frequency-based nonstationary dual wavelet frames in the 
distribution space @'(M), if and only if, 

OO 

JJ^W) + E4 fc (M) = 0, a.e. £eR,ke A\{0} (4.6) 
j=J 

(All the above infinite sums are in fact finite, since linij^ +00 Xjk = for all k G W.) and 

j'-i 

lim (^#(Aj-)+> J^j (Aj-) ) = 1 in the sense of distributions, (4.7) 

J'^+oo V ^ — 4 / 
3=J 

where A := U°° = j[Xj l Z] and 

r 

4(0:=E^^ + 27rA; )' fcGZ and 4(0 = 0, fcGE\Z (4.8) 
e=i 

and 

s j 

/£.(£) : = Y^^mW^i + 2 ^), keZ and 4i(0 = °> fc e ( 4 - 9 ) 
Proof. Let f , g G f^(R). By Lemma [3j we have 

EE( f '<;0, fc )(^ ; o, fc) g) = 2 7 r / E f(0g(£ + 2^)7^0 



(4.10) 



Since f and g have compact support and lim^ +00 Xj = 0, we observe that there exists an integer 

L-l' 



J fig such that f (f)g(f + 2ttA;) = for all k G [At 1 Z]\{0} and j ^ J f)g . Therefore, 

EE^ f -<:n>)(^:0.A.-g) = 2VT / f (0g(0^ (M) « V J > J fjg . (4.11) 

1=1 k&L ^ R 
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Sufficiency. For J' > J, define 

(f, g) := E E< f ' ^W<*l;0, fe , §) + E E E( f ' ''I't^mU,, S). 

1=1 k& j=J 1=1 keZ 

Therefore, by (|4.10|) . for J' > J, we have 

Sj' (f,g) = 2tt / E f (0g¥+2^)[4 Jfe (Aj0 + E4 fe (^0] (4-12) 

Now by (14.61) . for all J' > max(J, Jf ig ), we deduce that 

» J'-i 
Sj'(f,g) = 271 / f(Oi(Ok(AjO + E ^iOl (4-13) 

i=J 

Now by g2D and (Q3]l . we conclude that limj/_ H . 00 5^'(f, g) = 2tt j R f(f)I(0 = 27r ( f ; g>- 

Necessity. The proof of the necessity part is essentially the same as that of Lemma [5j Since 
linij^ +00 Xj = 0, the set A is discrete and closed. For any temporarily fixed (6l and k G A\{0}, 
it is important to notice that dist(/co, A\{fco}) > 0. Now the same argument as in the proof of 
Lemma [5] leads to (14.61) . Similarly, for any temporarily fixed (Gl, since dist(0, A\{0}) > 0, by 
(I4.13P , the same argument as in the proof of Lemma [5] leads to (14. 7p . □ 

For the particular case Xj = d - - 7 , the following result is a direct consequence of Theorem [T2l 

Corollary 13. Let d be an integer with |d| > 1 and J be an integer. Let <fr, $ in (14.1 ft and *Sfi, i&i 
in (14.21) be subsets of L l 2 oc (R) for all integers j is J. Then the pair in (14. 3p . with Xj = d~ J for all 
j ^ J, forms a pair of frequency- based nonstationary dual wavelet frames in the distribution space 
), if and only if, for all jo 6 N U {0} and k G Z\[dZ] ; 



r 30 s j+j 

E^(d J0 O^V°(£ + 2?rA; )) + EE^'' 7 ''^ <)V J " IJ (& :: J (£ + 27rjfc )) = 0, (4.14) 

1=1 j=0 1=1 

for almost every £ G M, and the following limit holds in the sense of distributions: 

r J'-l Sj 

jS^o (E^(^)^( d " J -) + E E^ ; ' (d ^ /(d J - ] ) = L ( 4 - 15 ) 

i=i j=j i=i 

Proof. Since Xj = d~ j , it is evident that flUD is just (l4~15|) . Note that A = U£Lj[cP'Z] = d J Z by 
d G Z. For any k G A\{0} = [d J Z]\{0}, we can write k = d J+JO /c for a unique integer j ^ and 
a unique k Q G Z\[dZ]. Now it is easy to check that (14.61) is equivalent to (I4.14p . □ 

As another application of Theorem [T2], using a similar argument as in the proof of Corollary ITOl 
we have the following result on frequency-based nonstationary tight wavelet frames in L 2 ( 



Corollary 14. Let J be an integer and {Xj} < jL J be a sequence of nonzero real numbers such that 
linij_^ +00 Xj = 0. Let <fr in (14.11) and SE^ in (14.21) be subsets of distributions in for all integers 

j ^ J . Then the following statements are equivalent 

(1) FWSj($; {W}f = j) is a frequency-based nonstationary tight wavelet frame in L 2 (M.), that 
is, <£, V j C L 2 {R) for all j ^ J and 

r oo Sj 

E E K f > *>Wl 2 + E E E K f > <;o, fc )! 2 = 2vr||f ||1 2(R) , V f G L 2 (R); (4.16) 
£=i fcez j=j £=i fcez 

(2) the pair (FWSj(3>; {^} < jLj)-, FWSj(<fr; {^}'jLj)) forms a pair of frequency-based nonsta- 
tionary dual wavelet frames in the distribution space @'(M); 

(3) * J C Lf, oc (R) and (TO . (j47D fcoW urato $ := $ and # J := * J /or a// j ^ J. 
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(4) there exist <p l , . . . ,<f) r , if?' 1 , . . . , ip j > Sj G L 2 (M) for all j ^ J such that (p l = (f) 1 ,. . . ,ip r 
0^ i/ji- 1 = ^3,1 : . . . 5 = ijjhsj for all j ^ J, and 



OO »j 



EEi</'^*>i 2 +EEEi</>^i sJfc >r = ii/iiU). v /^w ( 4 - 17 ) 

For frequency-based nonstationary wavelets with multiresolution-like structure, we have 

Proposition 15. Let J be an integer and {Xj}°^ Jo be a sequence of nonzero real numbers such 
that linij_+ +OC) Xj = 0. For integers j ^ Jo, /e£ & , \£ J in ( 14.2ft and 

$i : = {y? 3,1 , . . . , $ J := . . . , 0^} (4.18) 

be subsets of L l ° c (R). Then 

(FWSj($ J ; {v^'}~ j), FWSj($ J ; {* J }°° =J )) (4.19) 

forms a pair of frequency- based nonstationary dual wavelet frames in the distribution space f^'(R) 
for every integer J ^ J , if and only if, 

+ 4 fc ( A iO = o.e. £ G R, fc G [XfZ] U [At + \Z], j ^ J (4.20) 

and 

lim \ (p^' i (\j-)(p : '' i (Xj-) — 1 in the sense of distributions, (4-21) 

£=1 

where 1^, fcGf, are defined in (14 ,9p and 
r j 

4,(0 := £)^(O£^(£ + 27rA 0> and 4,(0=0, G R\Z. (4.22) 

£=1 

Proof. By the same argument as in Theorem El we see that the pair in (14. 19j) forms a pair of 
frequency-based nonstationary dual wavelet frames in S>'(M) for all integers J ^ Jo, if and only if, 



E E^ f ' V>£o,*> (^;0,fc> §) + E E^ f ' ^;0,fc) (^;0,fc> §) 

r j+ i ^ • ) 



and 

.lim ^X)< f '^f;o.*><^;o,fc.8> = Mf.g>, f,ge^(R). (4.24) 

7— >+oo -" ' ' 

£=i fcez 

By Lemma El (fl~20]) is equivalent to (I4T23D . By Lemma H, (T4~2Tj) is equivalent to (fl~24j) . □ 

For a(£) = YlkLk! a(k)e~ lk ^ with a{ki)a{k 2 ) ^ 0, the degree of a is defined to be deg(a) : = 
max(|A;i|, l^l). We finish this paper by the following result which connects a nonstationary wavelet 
filter bank obtained via a generalized nonstationary oblique extension principle with a pair of 
frequency-based nonstationary dual wavelet frames in the distribution space. 

Theorem 16. Let {dj}j° =1 be a sequence of nonzero integers such that hm.j^ +00 YYn=i \d n \ = °°- 

Define Xq := 1 and Xj := frin=i^«) f or a ^ 3 e Let a J and aP , j G N, be 2ir-periodic 
trigonometric polynomials such that a J (0) = a J (0) = 1 for all j G N and 

oo oo 

c := E l A J'l de g( aJ )ll ai IUoo(T) < oo and C:=J2 \ X j\ deg(a J ') || a J || ioo(T) < oo. (4.25) 
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Define 

oo oo 

<^(0 := II * J+n ( VnA7 x and := J] ^(X^Xj^), f6R,j6 NU{0}. (4.26) 

n=l n=l 

Tnen a// (p* -1 , l P J ~ 1 ,j G N are elements in L l ^(R) satisfying 

(pi- 1 ^) = «^ ^'"'(diO = ^'(0^(0, V^gKjJ'gN. (4.27) 

Lei 6P>\ . . . , 0?>i-i, b^'' 1 , . . . , fr^'- 1 and fik 1 , . . . , 0?>J-i, t>\ . . . , b^- 1 mt/i r^, Sj_i G N ana 1 
j G N fe 27r -periodic measurable functions in L 2 oc (IR). Define 

<^'- M (0 := ^(0^' _1 (0 «™d ^'-^(0 := ^(0^' _1 (0: ieN,^=l,...,r i _ 1 , (4.28) 
^ ; ''id,0: b-'-'SO^iO and ^ '-'(d/): b^(£)^'(£), J € N, £ = 1, . . . , Si _i. (4.29) 

Tnen <fr J , <fr J m ( 14. 18ft and VE^, ^ m (I4.2f) are subsets o/L 2 oc (R) /or a// j ^ 0. Moreover, the pair 
in (I4.19P forms a pair of frequency-based nonstationary dual wavelet frames in @'(M.) for every 
integer J ^ 0, if and only if, for all j G N, 

Sj-l 

9^(d,0^(0^(0 + ^Wffi&*{£) = 6^(0, a.e.e G R, (4.30) 

Sj-l 

^(djO^m^^ + 2 -ff) + ^WffiV*{£ + 2 -ff) = 0, a.e. e G R, w G Z\[d,Z], (4.31) 

and 

lim Q 3+ (A,-) = 1 in the sense of distributions, (4.32) 

wnere 

£=1 

Proof. We first establish an inequality on the decay of ip 3 , which plays a critical role in this proof 
to show g2ID. By aJ'(O) = #(0) = 1 and Bernstein inequality ||[a?]'|| Loof r) ^ deg(a')|| [a J ']|| Loo(T) , 

|l-a J '(OI = /V]'(CK ^ll[aTI|L 00 (T)|el^le|deg(^)||a^|| WT) , (4.34) 
./o 

where [a- 7 ]' denotes the derivative of a J . Observe a simple inequality 

kl ^ e 11 ^" 1 VzGC. (4.35) 

In fact, denote z = re td for r ^ and 6 G R. It is easy to prove that e 2r ^ (r + l) 2 for r ^ 0. 
Consequently, e 2 ' 2 ' ^ (r + l) 2 ^ r 2 — 2r cos# + 1 = \z — 1| 2 , which leads to (14.351) . 

Now for any mi, ni2 G N U {+00} with m\ ^ 777,2, by (I4.34p and (14.351) . we deduce that 

m 2 

j \\ a j+n (X j+n ^ e E — Jl-a^^+nOI ^ e klE™J mi |Ai+n|deg(ai+")||a^+»|| Loo(T) _ 
n=mi 

Hence, by ( I4.25p . we have 

J| a j+n (A i+ „0 ^ e c|fl V£ G R,mi,m 2 G NU {+00} with m x ^ m 2 . (4.36) 

n=mi 

On the other hand, by a similar idea as in [2"2l Page 93] and |26, Page 932], we have 

7712 T,2 7T72 

1- a*«(A i+B fl= £ (l - a^(A J+m 0) ( a i+ "(A;+n£)) , 

n=mi 7Tl=?77l 71=777 + 1 
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where n™=m 2 +i := 1- By (14.341) and (14.361) . we deduce from the above identity that 

m2 TU2 1712 

1-fl a^(A i+ nO|^^l £ |l-a J+m (A J+m Ol^e^l|e| £ |A i+m | deg(a^ m )||a^ m || Loo(T) . 

n=mi m=m\ m=m\ 

That is, for all m 1; m 2 G N U {+00} with m\ ^ m 2 , we have 

r«2 j+m2 

1- H^ +n (X J+n O\^e c ^\ £ |A m |deg(a m )||a m || ioo(T) V^Gl. (4.37) 

n=mi m=j+mi 

The above inequality implies the uniform convergence of YY^=i a,J+ri (Ay+nO f° r £ 011 an y bounded 
set. Since <^'(AjO = Un=i a j+n (A j+ „0, we conclude from lOBjl that <p J ' G L l £(R). Since 
0J> f 5 b?^ G L 2 oc (M), it is evident that all ip 3 ^,^^ are elements in L 2 oc (R). Similarly, we can 
prove (pi G L^f(R) and all (p?£,i\p£ are elements in L 2 oc (R). 

Note that ^-'(Xj-iO = a'(A/)^(A/). that is, ^(O = a% '(V^ '0- So, g^Z) holds. 
Also note that lim^oo Aj = and = for at most countably many £ G R. Note that 

(I4.20p with j being replaced by j — 1 is equivalent to 

E^ rT7 id~0^' -'(d/ + 27rdjife + 2vrw) + ^ V' J ~ 1/ (djO^" 1/ (d^ + 2nd jk + 2ttcu) 

(4.38) 



d, 

1=1 



for all G Z and = 0, . . . , |d 3 -| — 1, where 5(0) = 1 and 5 (a;) = for all 10 7^ 0. Now it is easy to 
directly verify that fl4~30]) and fOTD are equivalent to (OOj) . 

We now show that (OB is equivalent to (Q2]l . By ^'(M) = Il~=i a i+n (A i+n £) and (Q7|l . 
we have 

00 

ll-^-OKe^l £ |A m |deg(a m )||a m || Loo(T) , ^1,j6NU{0}. (4.39) 

m=j+l 

By a similar argument, we also have 

00 

£ |A m |deg(a m )||a m || Loo(T) , £gR,jgNU{0}. (4.40) 

m=j+l 

By ( 14.251) . we infer that 



lim V |A m |deg(a m )||a m || ioo(T) =0 and lim V |A m | deg(a m ) ||a m || Loo(T) = 0. 

j— >+oo ' * ' J 

m=j+l m=j'+l 

By (BHD , we deduce that 

E^wy^(A,o = e^^o^wy^o. 

£=1 

Now by (I4.39P and (14.401) . using Lebesgue dominated convergence theorem, by the same argument 
as in Theorem [2} we conclude that (14.211) holds if and only if (14.321) holds. By Proposition [151 the 
proof is completed. □ 
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